<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<title>nih: sampling/halton.h Source File</title>

<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
  $(document).ready(initResizable);
</script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
  $(document).ready(function() { searchBox.OnSelectItem(0); });
</script>

</head>
<body>
<div id="top"><!-- do not remove this div! -->


<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  
  
  <td style="padding-left: 0.5em;">
   <div id="projectname">nih
   &#160;<span id="projectnumber">0.5</span>
   </div>
   
  </td>
  
  
  
 </tr>
 </tbody>
</table>
</div>

<!-- Generated by Doxygen 1.7.5.1 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="index.html"><span>Main&#160;Page</span></a></li>
      <li><a href="modules.html"><span>Modules</span></a></li>
      <li><a href="annotated.html"><span>Classes</span></a></li>
      <li class="current"><a href="files.html"><span>Files</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
        <span class="left">
          <img id="MSearchSelect" src="search/mag_sel.png"
               onmouseover="return searchBox.OnSearchSelectShow()"
               onmouseout="return searchBox.OnSearchSelectHide()"
               alt=""/>
          <input type="text" id="MSearchField" value="Search" accesskey="S"
               onfocus="searchBox.OnSearchFieldFocus(true)" 
               onblur="searchBox.OnSearchFieldFocus(false)" 
               onkeyup="searchBox.OnSearchFieldChange(event)"/>
          </span><span class="right">
            <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
          </span>
        </div>
      </li>
    </ul>
  </div>
  <div id="navrow2" class="tabs2">
    <ul class="tablist">
      <li><a href="files.html"><span>File&#160;List</span></a></li>
    </ul>
  </div>
</div>
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
  initNavTree('halton_8h.html','');
</script>
<div id="doc-content">
<div class="header">
  <div class="headertitle">
<div class="title">sampling/halton.h</div>  </div>
</div>
<div class="contents">
<div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="preprocessor">#pragma once</span>
<a name="l00002"></a>00002 <span class="preprocessor"></span>
<a name="l00003"></a>00003 <span class="preprocessor">#include &lt;<a class="code" href="vector_8h.html" title="Define linear-algebra vector classes.">nih/linalg/vector.h</a>&gt;</span>
<a name="l00004"></a>00004 <span class="preprocessor">#ifdef __CUDACC__</span>
<a name="l00005"></a>00005 <span class="preprocessor"></span><span class="preprocessor">#include &lt;cuda_runtime_api.h&gt;</span>
<a name="l00006"></a>00006 <span class="preprocessor">#endif</span>
<a name="l00007"></a>00007 <span class="preprocessor"></span>
<a name="l00008"></a>00008 <span class="keyword">namespace </span>nih {
<a name="l00009"></a>00009 
<a name="l00010"></a>00010 <span class="keyword">namespace </span>halton {
<a name="l00011"></a>00011 <span class="keyword">extern</span> uint32 s_halton_bases[2000];
<a name="l00012"></a>00012 <span class="keyword">extern</span> <span class="keywordtype">float</span>  s_halton_sigma[];
<a name="l00013"></a>00013 } <span class="comment">// namespace halton</span>
<a name="l00014"></a>00014 
<a name="l00015"></a>00015 <span class="keyword">inline</span> NIH_HOST_DEVICE uint32 halton_bases(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> host_space tag)   { <span class="keywordflow">return</span> halton::s_halton_bases[b]; }
<a name="l00016"></a>00016 <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> host_space tag)   { <span class="keywordflow">return</span> halton::s_halton_sigma[b]; }
<a name="l00017"></a>00017 <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> int32 bj, <span class="keyword">const</span> host_space tag)     { <span class="keywordflow">return</span> halton::s_halton_sigma[b] / float(bj); }
<a name="l00018"></a>00018 
<a name="l00019"></a>00019 <span class="preprocessor">#ifdef __CUDACC__</span>
<a name="l00020"></a>00020 <span class="preprocessor"></span>
<a name="l00021"></a>00021 texture&lt;float, 1&gt;   tex_halton_sigma;
<a name="l00022"></a>00022 __constant__ uint32 c_halton_bases[2000];
<a name="l00023"></a>00023 
<a name="l00024"></a>00024 <span class="preprocessor">#if __CUDA_ARCH__</span>
<a name="l00025"></a>00025 <span class="preprocessor"></span>
<a name="l00026"></a>00026 <span class="keyword">inline</span> __host__ __device__ uint32 halton_bases(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> c_halton_bases[b]; }
<a name="l00027"></a>00027 <span class="keyword">inline</span> __host__ __device__ <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> tex1Dfetch( tex_halton_sigma, b ); }
<a name="l00028"></a>00028 <span class="keyword">inline</span> __host__ __device__ <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> int32 bj, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> __fdividef( tex1Dfetch( tex_halton_sigma, b ), <span class="keywordtype">float</span>(bj) ); }
<a name="l00029"></a>00029 
<a name="l00030"></a>00030 <span class="preprocessor">#else</span>
<a name="l00031"></a>00031 <span class="preprocessor"></span>
<a name="l00032"></a>00032 <span class="keyword">inline</span> __host__ __device__ uint32 halton_bases(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> 0; }
<a name="l00033"></a>00033 <span class="keyword">inline</span> __host__ __device__ <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> 0.0f; }
<a name="l00034"></a>00034 <span class="keyword">inline</span> __host__ __device__ <span class="keywordtype">float</span>  halton_sigma(<span class="keyword">const</span> int32 b, <span class="keyword">const</span> int32 bj, <span class="keyword">const</span> device_space tag) { <span class="keywordflow">return</span> 0.0f; }
<a name="l00035"></a>00035 
<a name="l00036"></a>00036 <span class="preprocessor">#endif</span>
<a name="l00037"></a>00037 <span class="preprocessor"></span>
<a name="l00038"></a>00038 <span class="preprocessor">#endif</span>
<a name="l00039"></a>00039 <span class="preprocessor"></span>
<a name="l00043"></a>00043 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag = host_space&gt;
<a name="l00044"></a><a class="code" href="classnih_1_1_halton__sampler.html">00044</a> <span class="keyword">class </span><a class="code" href="classnih_1_1_halton__sampler.html">Halton_sampler</a>
<a name="l00045"></a>00045 {
<a name="l00046"></a>00046 <span class="keyword">public</span>:
<a name="l00047"></a>00047         <span class="keyword">inline</span> NIH_HOST_DEVICE <a class="code" href="classnih_1_1_halton__sampler.html">Halton_sampler</a> ();
<a name="l00048"></a>00048           <span class="comment">// Constructs an empty sampling context.</span>
<a name="l00049"></a>00049           <span class="comment">// Child contexts will inherit the randomization seed.</span>
<a name="l00050"></a>00050 
<a name="l00051"></a>00051         <span class="keyword">inline</span> NIH_HOST_DEVICE ~Halton_sampler ();
<a name="l00052"></a>00052           <span class="comment">// Destructor</span>
<a name="l00053"></a>00053 
<a name="l00054"></a>00054         <span class="keyword">inline</span> NIH_HOST_DEVICE Halton_sampler split (<span class="keyword">const</span> uint32 dim, <span class="keyword">const</span> uint32 n, <span class="keyword">const</span> uint32 start = 0);
<a name="l00055"></a>00055           <span class="comment">// Dependent splitting function, creating a child sampling</span>
<a name="l00056"></a>00056           <span class="comment">// context with &#39;dim&#39; dimensions and &#39;n&#39; samples.</span>
<a name="l00057"></a>00057 
<a name="l00058"></a>00058     <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keywordtype">bool</span>  next_sample(<span class="keywordtype">float</span>* u);
<a name="l00059"></a>00059     <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keywordtype">float</span> sample(<span class="keyword">const</span> int32 i);
<a name="l00060"></a>00060 
<a name="l00061"></a>00061     <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keywordtype">bool</span>  advance() { <span class="keywordflow">return</span> ++dLI &lt; dLN; }
<a name="l00062"></a>00062           <span class="comment">// Next sample from the sampling loop. Returns 0 if the</span>
<a name="l00063"></a>00063           <span class="comment">// loop is already finished, 1 otherwise.</span>
<a name="l00064"></a>00064 
<a name="l00065"></a>00065         <span class="keyword">inline</span> <span class="keyword">static</span> <span class="keywordtype">void</span> init(<span class="keywordtype">void</span>);
<a name="l00066"></a>00066           <span class="comment">// Initialization function</span>
<a name="l00067"></a>00067 
<a name="l00068"></a>00068         <span class="keyword">inline</span> NIH_HOST_DEVICE <span class="keyword">static</span> <span class="keywordtype">float</span> phi (uint32 n, uint32 b);
<a name="l00069"></a>00069           <span class="comment">// Radical inverse function</span>
<a name="l00070"></a>00070 
<a name="l00071"></a>00071         <span class="keyword">static</span> <span class="keyword">const</span> uint32  sMAXDIM = 200;
<a name="l00072"></a>00072     <span class="keyword">static</span> <span class="keyword">const</span> uint32  sBASESCOUNT = 9974;
<a name="l00073"></a>00073 
<a name="l00074"></a>00074 <span class="keyword">private</span>:
<a name="l00075"></a>00075         uint32 dGI; <span class="comment">// Global instance number</span>
<a name="l00076"></a>00076         uint32 dGD; <span class="comment">// Global integral dimension</span>
<a name="l00077"></a>00077         uint32 dLI; <span class="comment">// Local instance number</span>
<a name="l00078"></a>00078         uint32 dLD; <span class="comment">// Local integral dimension</span>
<a name="l00079"></a>00079         uint32 dLN; <span class="comment">// Local number of samples</span>
<a name="l00080"></a>00080 };
<a name="l00081"></a>00081 
<a name="l00082"></a>00082 <span class="preprocessor">#ifdef __CUDACC__</span>
<a name="l00083"></a>00083 <span class="preprocessor"></span><span class="keyword">struct </span>Halton_sampler_binder
<a name="l00084"></a>00084 {
<a name="l00085"></a>00085     Halton_sampler_binder()
<a name="l00086"></a>00086     {
<a name="l00087"></a>00087         <span class="comment">// copy bases &amp; sigma to constant memory</span>
<a name="l00088"></a>00088         cudaMemcpyToSymbol( c_halton_bases, halton::s_halton_bases, <span class="keyword">sizeof</span>(uint32)*2000 );
<a name="l00089"></a>00089 
<a name="l00090"></a>00090         <span class="keyword">const</span> uint32 bases_count = <a class="code" href="classnih_1_1_halton__sampler.html">Halton_sampler&lt;device_space&gt;::sBASESCOUNT</a>;
<a name="l00091"></a>00091 
<a name="l00092"></a>00092         cudaMalloc( (<span class="keywordtype">void</span>**)&amp;g_sigma, <span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)*bases_count*bases_count );
<a name="l00093"></a>00093         <span class="comment">//cudaMalloc( (void**)&amp;g_bases, sizeof(uint32)*2000 );</span>
<a name="l00094"></a>00094         <span class="keywordflow">if</span> (g_sigma == NULL <span class="comment">/*|| g_bases == NULL*/</span>)
<a name="l00095"></a>00095         {
<a name="l00096"></a>00096             fprintf(stderr, <span class="stringliteral">&quot;Halton_sampler_binder(): allocation failed\n&quot;</span>);
<a name="l00097"></a>00097             <span class="keywordflow">return</span>;
<a name="l00098"></a>00098         }
<a name="l00099"></a>00099         cudaMemcpy( g_sigma, halton::s_halton_sigma, <span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)*bases_count*bases_count, cudaMemcpyHostToDevice );
<a name="l00100"></a>00100         <span class="comment">//cudaMemcpy( g_bases, halton::s_halton_bases, sizeof(int32)*2000, cudaMemcpyHostToDevice );</span>
<a name="l00101"></a>00101 
<a name="l00102"></a>00102         cudaChannelFormatDesc channel_desc = cudaCreateChannelDesc&lt;float&gt;();
<a name="l00103"></a>00103         cudaBindTexture( 0, &amp;tex_halton_sigma, g_sigma, &amp;channel_desc, bases_count*bases_count );
<a name="l00104"></a>00104     }
<a name="l00105"></a>00105     ~Halton_sampler_binder()
<a name="l00106"></a>00106     {
<a name="l00107"></a>00107         <span class="comment">//cudaUnbindTexture( &amp;tex_halton_sigma );</span>
<a name="l00108"></a>00108         cudaFree( g_sigma );
<a name="l00109"></a>00109         cudaFree( g_bases );
<a name="l00110"></a>00110     }
<a name="l00111"></a>00111 
<a name="l00112"></a>00112     <span class="keywordtype">float</span>* g_sigma;
<a name="l00113"></a>00113     int32* g_bases;
<a name="l00114"></a>00114 };
<a name="l00115"></a>00115 <span class="preprocessor">#endif</span>
<a name="l00116"></a>00116 <span class="preprocessor"></span>
<a name="l00117"></a>00117 <span class="comment">// ============================================================================ //</span>
<a name="l00118"></a>00118 <span class="comment">//                             IMPLEMENTATION                                   //</span>
<a name="l00119"></a>00119 <span class="comment">// ============================================================================ //</span>
<a name="l00120"></a>00120 
<a name="l00121"></a>00121 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00122"></a>00122 <span class="keyword">inline</span> Halton_sampler&lt;space_tag&gt;::Halton_sampler () :
<a name="l00123"></a>00123         dGI ( 0 ), dGD( 0 ), dLI( 0 ), dLD( 0 ), dLN( 0 )
<a name="l00124"></a>00124 {
<a name="l00125"></a>00125 }
<a name="l00126"></a>00126 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00127"></a>00127 <span class="keyword">inline</span> Halton_sampler&lt;space_tag&gt;::~Halton_sampler ()
<a name="l00128"></a>00128 {
<a name="l00129"></a>00129 }
<a name="l00130"></a>00130 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00131"></a>00131 <span class="keyword">inline</span> Halton_sampler&lt;space_tag&gt; Halton_sampler&lt;space_tag&gt;::split (<span class="keyword">const</span> uint32 dim, <span class="keyword">const</span> uint32 n, <span class="keyword">const</span> uint32 start)
<a name="l00132"></a>00132 {
<a name="l00133"></a>00133         Halton_sampler ctxt;
<a name="l00134"></a>00134         ctxt.dGI = dGI + dLI; <span class="comment">// Decorrelate dependent samples</span>
<a name="l00135"></a>00135         ctxt.dGD = dGD + dim; <span class="comment">// Compute global dimension</span>
<a name="l00136"></a>00136         ctxt.dLI = start;
<a name="l00137"></a>00137         ctxt.dLD = dim;
<a name="l00138"></a>00138         ctxt.dLN = n + start;
<a name="l00139"></a>00139         <span class="keywordflow">return</span> ctxt;
<a name="l00140"></a>00140 }
<a name="l00141"></a>00141 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00142"></a>00142 <span class="keyword">inline</span> <span class="keywordtype">bool</span> Halton_sampler&lt;space_tag&gt;::next_sample(<span class="keywordtype">float</span>* u)
<a name="l00143"></a>00143 {
<a name="l00144"></a>00144         <span class="keywordflow">if</span> (dLI == dLN) <span class="comment">// Check if loop is finished</span>
<a name="l00145"></a>00145                 <span class="keywordflow">return</span> <span class="keyword">false</span>;
<a name="l00146"></a>00146 
<a name="l00147"></a>00147     <span class="keywordflow">for</span> (uint32 j = 0; j &lt; dLD; ++j)
<a name="l00148"></a>00148         u[j] = sample( j );
<a name="l00149"></a>00149 
<a name="l00150"></a>00150     <span class="keywordflow">return</span> advance();
<a name="l00151"></a>00151 }
<a name="l00152"></a>00152 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00153"></a>00153 <span class="keyword">inline</span> <span class="keywordtype">float</span> Halton_sampler&lt;space_tag&gt;::sample (<span class="keyword">const</span> int32 j)
<a name="l00154"></a>00154 {
<a name="l00155"></a>00155         <span class="comment">// Build sample point</span>
<a name="l00156"></a>00156     <span class="keyword">const</span> <span class="keywordtype">float</span> r = dGI ? 
<a name="l00157"></a>00157         fmodf(
<a name="l00158"></a>00158             phi( dGI, dGD - dLD + j ) +
<a name="l00159"></a>00159             phi( dLI, j ), 1.0f ) :
<a name="l00160"></a>00160             phi( dLI, j );
<a name="l00161"></a>00161     <span class="keywordflow">return</span> r;
<a name="l00162"></a>00162 }
<a name="l00163"></a>00163 
<a name="l00164"></a>00164 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00165"></a>00165 <span class="keyword">inline</span> <span class="keywordtype">float</span> Halton_sampler&lt;space_tag&gt;::phi (uint32 n, uint32 bindex)
<a name="l00166"></a>00166 {
<a name="l00167"></a>00167         <span class="keywordtype">float</span>  result = 0.0;
<a name="l00168"></a>00168         uint32  remainder;
<a name="l00169"></a>00169         uint32  m, bj = 1, b = halton_bases( bindex, space_tag() );
<a name="l00170"></a>00170 
<a name="l00171"></a>00171         <span class="keywordflow">do</span>
<a name="l00172"></a>00172         {
<a name="l00173"></a>00173                 bj *= b;
<a name="l00174"></a>00174                 m   = n;
<a name="l00175"></a>00175                 n  /= b;
<a name="l00176"></a>00176 
<a name="l00177"></a>00177                 remainder = m - n * b;
<a name="l00178"></a>00178 
<a name="l00179"></a>00179                 result += halton_sigma( b*sBASESCOUNT + remainder, bj, space_tag() );
<a name="l00180"></a>00180         } <span class="keywordflow">while</span> (n &gt; 0);
<a name="l00181"></a>00181 
<a name="l00182"></a>00182         <span class="keywordflow">return</span> result;
<a name="l00183"></a>00183 }
<a name="l00184"></a>00184 
<a name="l00185"></a>00185 <span class="preprocessor">#define SIGMA(b,i) halton::s_halton_sigma[ ((b) * sBASESCOUNT + (i)) ]</span>
<a name="l00186"></a>00186 <span class="preprocessor"></span>
<a name="l00187"></a>00187 <span class="keyword">template</span> &lt;<span class="keyword">typename</span> space_tag&gt;
<a name="l00188"></a>00188 <span class="keyword">inline</span> <span class="keywordtype">void</span> Halton_sampler&lt;space_tag&gt;::init (<span class="keywordtype">void</span>)
<a name="l00189"></a>00189 {
<a name="l00190"></a>00190         <span class="comment">// Build Faure permutation tables</span>
<a name="l00191"></a>00191         SIGMA(2,0) = 0;
<a name="l00192"></a>00192         SIGMA(2,1) = 1;
<a name="l00193"></a>00193         <span class="keywordflow">for</span> (uint32 b = 3; b &lt; sBASESCOUNT; b++)
<a name="l00194"></a>00194         {
<a name="l00195"></a>00195                 uint32 i;
<a name="l00196"></a>00196 
<a name="l00197"></a>00197                 <span class="keywordflow">if</span> (b % 2 == 0)
<a name="l00198"></a>00198                 {
<a name="l00199"></a>00199                         uint32 bByHalf = b / 2;
<a name="l00200"></a>00200 
<a name="l00201"></a>00201                         <span class="keywordflow">for</span> (i = 0; i &lt; bByHalf; i++)
<a name="l00202"></a>00202                                 SIGMA( b, i ) = 2 * SIGMA( bByHalf, i );
<a name="l00203"></a>00203 
<a name="l00204"></a>00204                         <span class="keywordflow">for</span> (i = 0; i &lt; bByHalf; i++)
<a name="l00205"></a>00205                                 SIGMA( b, bByHalf + i ) = 2 * SIGMA( bByHalf, i ) + 1;
<a name="l00206"></a>00206                 }
<a name="l00207"></a>00207                 <span class="keywordflow">else</span>
<a name="l00208"></a>00208                 {
<a name="l00209"></a>00209                         uint32 bMinusOneByHalf = (b - 1) / 2;
<a name="l00210"></a>00210 
<a name="l00211"></a>00211                         <span class="keywordflow">for</span> (i = 0; i &lt; bMinusOneByHalf; i++)
<a name="l00212"></a>00212                         {
<a name="l00213"></a>00213                                 <span class="keywordtype">float</span> s = SIGMA( b-1, i );
<a name="l00214"></a>00214                                 SIGMA( b, i ) = s &lt; bMinusOneByHalf ? s : s+1;
<a name="l00215"></a>00215                         }
<a name="l00216"></a>00216 
<a name="l00217"></a>00217                         <span class="keywordflow">for</span> (i = bMinusOneByHalf; i &lt; b-1; i++)
<a name="l00218"></a>00218                         {
<a name="l00219"></a>00219                                 <span class="keywordtype">float</span> s = SIGMA( b-1, i );
<a name="l00220"></a>00220                                 SIGMA( b, i+1 ) = s &lt; bMinusOneByHalf ? s : s+1;
<a name="l00221"></a>00221                         }
<a name="l00222"></a>00222 
<a name="l00223"></a>00223                         SIGMA( b, bMinusOneByHalf ) = (float)bMinusOneByHalf;
<a name="l00224"></a>00224                 }
<a name="l00225"></a>00225         }
<a name="l00226"></a>00226 }
<a name="l00227"></a>00227 
<a name="l00228"></a>00228 <span class="preprocessor">#undef SIGMA</span>
<a name="l00229"></a>00229 <span class="preprocessor"></span>
<a name="l00230"></a>00230 } <span class="comment">// namespace nih</span>
</pre></div></div>
</div>
  <div id="nav-path" class="navpath">
    <ul>
      <li class="navelem"><b>halton.h</b>      </li>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
     onmouseover="return searchBox.OnSearchSelectShow()"
     onmouseout="return searchBox.OnSearchSelectHide()"
     onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark">&#160;</span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark">&#160;</span>Classes</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark">&#160;</span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark">&#160;</span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark">&#160;</span>Variables</a></div>

<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0" 
        name="MSearchResults" id="MSearchResults">
</iframe>
</div>


    <li class="footer">Generated on Thu Dec 15 2011 11:56:20 for nih by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.5.1 </li>
   </ul>
 </div>


</body>
</html>
